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Two idealized models for the preheat stage of liquid droplets are analyzed theoretically. These models contain 
the effects of transient heat conduction and evaporation. It is assumed that the droplet surface area decreases 
linearly with time. This assumption necessitates the solution of moving boundary problems. These models, 
however, do not consider gas-phase mass transport. In the finite-gradient model, the temperatures of both the 
droplet and surrounding hot gases vary spatially and temporally. In the zero-gradient model the gas temperature 
varies spatially and temporally but the droplet temperature varies only temporally, i.e., the droplet temperature is 
spatially uniform. Numerical examples, which require extensive calculations of confluent hypergeometric func- 
tions, are presented for typical values of the droplet latent heat and evaporation rate constant. The temperature 
profiles given by the finite-gradient and zero-gradient models agree to within 20 percent of each other for all cases 
examined. 
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1 . Introduction 

The ignition of a droplet of conventional fuel consists 
essentially of two stages [I] 1 . During the first or preheat stage 
heat flows from the hot surrounding gas to the droplet, 
causing the droplet temperature to rise and liquid fuel to 
evaporate from the droplet surface. In the second stage, 
ignition occurs in the gaseous mixture of fuel and oxidizer 
surrounding the droplet. The preheat stage is dominated by 
transient processes and provides the motivation for the sub- 
ject of this paper. 

Wise and Ablow [2], Parks, et al. [3], and Waldman [4] 
have previously analyzed the effects of transient droplet 
heating by neglecting the existence of internal circulation and 
by neglecting selected terms in the heat conduction equa- 
tions. In addition they assume that the droplet surface re- 
gresses linearly with time. In a previous publication, Bennett 
and Kayser [5] included the effects of internal circulation but 
also neglected evaporation and regression of the droplet 
surface. We now include these effects in the present paper. 

The neglect of internal circulation can in many instances 
lead to gross underestimates of the heat transfer rates within 
the droplet. Available experimental evidence indicates that 
in most cases the heat transfer rate inside the droplet is much 
faster than is possible with heat conduction alone. El Wakil, 
et al. [6] observed vigorous circulation within the droplet and 
showed that the droplet temperature is uniform even during 
the initial transient period. This indicates that the rate at 
which heat is transferred to the droplet is much slower than 
the rate of internal mixing. As pointed out by Law [7], 
assuming the droplet temperature to be uniform in a sense 
circumvents the difficult task of describing internal circula- 
tion but still includes the effects due to its presence; that is, 
uniform droplet temperature represents the limit of rapid 
internal circulation. 

Two models for the transient preheat stage are presented. 
In the first model, called the finite-gradient model, the 
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effects of transient heat conduction and evaporation with 
linear surface regression rate are analyzed. The droplet and 
gas temperatures vary spatially and temporally and the effects 
of internal circulation are neglected. In the second model, 
called the zero-gradient model, we include the effects of 
transient heat flow, evaporation, and internal circulation. 
The latter is included by assuming the droplet temperature to 
be temporally varying but spatially uniform as discussed 
above. The gas temperature is again both time and space 
dependent. 

The results of the finite-gradient and zero-gradient models 
are applied to the preheat stage of droplets in hot gases. The 
temperature profiles are used to predict the time duration 
between introduction of the droplet into a hot atmosphere and 
its reaching a given temperature. In addition we compare the 
analogous times predicted by a particularly simple and accu- 
rate approximation to the exact expression derived for the 
zero-gradient model. 

In the next section we describe the physical assumptions 
for the preheat stage common to both models. Section 3 
contains the mathematical description of these assumptions. 
We derive in section 4 the temperature profiles of the droplet 
and gas predicted by the finite-gradient model while in 
section 5 we give the corresponding expressions for the zero- 
gradient model. We present in the last section numerical 
examples for both models. The temperature profiles obtained 
are very similar for both models and are within 20 percent of 
each other in all cases examined. 

2. Preheat Stage 

Williams has examined the problem of the assumptions 
underlying droplet vaporization and combustion theories in 
some detail [8]. In models for preheating, many researchers 
consider the droplet to be a sphere with temperature inde- 
pendent thermal properties and neglect diffusion effects [1]. 
Others include thermal diffusion effects but later in their 
theoretical analyses neglect some terms in the diffusion equa- 
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tions [3]. In addition, most theoretical investigators neglect 
the influence of gravity, forced convection, and radiant heat 
transfer. Combustion chambers and experimental techniques 
such as freely falling droplets exist to eliminate gravitational 
effects. Proper selection of some experimental parameters 
such as ambient pressure and initial droplet size may mini- 
mize the other effects not included in the model for preheat- 
ing. Several researchers have shown that except for the very 
heavy fuel oils, radiant heat transfer from the hot gas or from 
adjacent droplets is negligible [7, 9]. 

The models examined in this paper contain only thermal 
diffusion effects. We assume that a spherical liquid droplet of 
initial radius r = r d0 is inserted at time t = into a hot 
atmosphere (e.g., air). The initial temperature of the droplet 
at time t = is T d0 - The hot gas is unbounded and initially at 
a constant temperature T qq. For t > the temperature at 
infinity is kept at Tg . 

The inclusion of evaporation demands that the droplet 
radius change with time. In the present work a linear surface 
regression rate is assumed: 



R 2 (t) = 4) = Kt 



(i) 



R(t) is a direct specification of the motion of the interface r = 
R(t) with time and k e denotes the evaporation rate constant. A 
self-contained theory would obtain R(t) indirectly by adding 
to the boundary conditions at the interface the requirement 
that the local vapor pressure be determined in the solution 
process. A moving boundary problem in which R(t) is directly 
specified falls into the general class of problems called 
inverse-Stefan problems while those in which it is determined 
in the solution process are called Stefan problems [10]. It 
must be stressed at this point that the solutions which we will 
obtain will be applicable only in those time regimes for which 
eq (1) is an accurate representation of the actual surface 
regression rate. 

We denote the density, specific heat, and the thermal 
conductivity of the droplet respectively by d d , C d , andK^ and 
of the gas respectively by dg, C g and K g . We assume that 
these thermal properties remain spatially and temporally 
constant during preheat. In addition, the ambient pressure of 
the host gas does not change. Using these and the above 
assumptions we seek to compute the dependence of the 
temperature T(r, t) upon the radial distance r and time t. In 
particular, we shall compute the time required for a droplet to 
reach a given temperature. 

3. Theoretical Analysis 

We first consider the finite-gradient model. The distribu- 
tion of temperature within the system is governed by the 
Fourier heat conduction equation for the case in which the 
isothermal surfaces are concentric spheres. Let T d (r, t), 
T g (r, t), a d = (KJddCd), and a g = {K g /dgC g ) denote the 
temperatures and thermal diffusivities of the droplet and 
surrounding gas respectively. The governing equation fort > 
and < r < R (inside the droplet) reads 



dT d (r, t) = ^/ d 2 T d (r,t) + 2 &T jr, t) * 



dt 



dr 2 



dr 



and for t ^ and r > R (outside the droplet) 



dT g (r, t) 
dt 



d 2 T g (r, t) 2 3T£r, t) 



dr 2 



+ 



dr 



(3) 



The quantity R(t) specifies the motion of the interface with 
time and is given in eq. (1). 

The specification of the boundary conditions completes the 
statement of the problem. The temperature T(r, t) has the 
form 



T d (r, t) = T d0 for < r < R and t < 



and 



T g (r, t) = T g0 for r > R and t 
and it becomes for t > 



(4) 



(5) 



7V, t) 



T d (r, t) for < r < R 
T g (r, t) for r > R. 



The temperature is finite everywhere, namely 7^(0, t) is finite 
and 



lim T,(r, t) = % 



go 



(6) 



The symmetry condition that no heat flux exists at the center 
of the droplet is 



dT d (0, t) 
'dr 



= 



(7) 



dT d (0, t)/dr is of course taken to mean the partial derivative 
of T d (r, t) which respect to r evaluated at r — 0. The problem 
statement for the finite-gradient model is completed by two 
conditions at the interface r = R which match the tempera- 
ture and heat flux 



T d (R~, t) = T g (R + , t) 



K, 



t>T d (R-, t) 



K. 



dT g {R + , t) 



+ Ldda 



dR(t) 

dt 



(8) 

(9) 



where R ± = lim (R ± Q and L d is the latent heat of vapori- 

zation per unit mass of the condensed material. From eq (1) 
we note that dR/dt goes off to infinity near the end of the 
droplet lifetime so that the difference between the thermal 
conductivities times the temperature gradients must become 
unbounded. More will be said about this point in subsequent 
sections. 

We find it convenient at this point to introduce dimension- 
less quantities which in this paper will be denoted by Greek 
letters. They are defined as follows: 

(1) The time T is chosen to move the singular point of 
complete droplet evaporation to infinity, r = 2{ajk e ) In 
(r d0 /R(t)). As the evaporation rate constant k e goes to zero r 

(2) reduces to T = a d t/r d{) . 

(2) The distance coordinate f] is chosen to remove the time 
dependence of the position of the boundary, 7] = r/R(t). 
Distance is thus measured in instantaneous droplet radii. 



258 



When k € = 0, T) = i)o = r/r d0 . 

(3). The dimensionsless temperature ratios are defined by 

8/77, r) = (Tj(r, /) - T do )/(Tgo - F d0 ), J = dorg. 

6 rf (7j, r) is the fraction of the maximum possible temperature 
rise. In addition, we define a = a d /a gi /3 = A^/A,,. e = 
k e /a%, and k = L d /\C d (Tgo - Too)]. 



With the above definitions, eqs (2)-(9) become 



ae d (T), t) = d^n, r) 

dr drj 2 



+ 



(- 



2 6T7\ ae d (T7, t) 
2 / a?? 



(0 < 77 < 1) (10) 



a 



^(77, r) _ 3*6,(7), r) 



dr 



+ 



'2 a 2 erA f)B,( 77 , t) 



/2 areTAc 
U " ~~2/ 



dr) 

e d (77, o) = o 

6,(77, 0) - 1 
lim 60(17, r) = 1 

7)—>ao 



(7) > 1) (ID 

(12) 
(13) 

(14) 



must equal the heat flux from the gas. Equations (2) and (9) 
are thus replaced by the following equality 



K a r L d d (l 



3 



f// 



dr 



d* 



(19) 



Equations (3)-(8) with T d (r, t) and 7^(r, /) replaced respec- 
tively by T dz (t) and 7 gz (r, t) and eq (19) are the mathematical 
statement of the zero-gradient model. 

In terms of the dimensionless quantities introduced ear- 
lier, the working equations in the zero-gradient model are 



. 2 de B :( V , t) d 2 OUr), t) 



3t 



+ 



drf 



2 a 2 er,\;)Q az (r), 7) 



/2 a^erA 
U~ ~~2/ 



dr) 



6^(0) = 



e„(ij. 0) = 1 



lim 692(17, t) = 1 



e rf2 (T) = e 9Z (i + , r) 



(tj>1) (20) 



(21) 



(22) 



(23) 



(24) 



ae d (o, T ) 
e d (i- r)= e fl (i + , T ) 



^i-ri^ae^.T),^ 



cJtj 



dT} 



(15) 



(16) 



(17) 



3 rfr 



^e g ,(i + , t) 

C*77 



2 



(25) 



In the next two sections we solve the equations for the 
finite-gradient and zero-gradient models by taking their La- 
place transforms. The Laplace transform of the reduced 
temperature 6(77, r) is denoted by 



Equations (10)-(17) are the complete statement in dimen- 
sionless form for the finite-gradient model. 

We now consider the zero-gradient model. As was men- 
tioned in section 1, the internal circulation inside the droplet 
may be sufficient in many cases to maintain a spatially 
uniform temperature within the droplet. When such condi- 
tions prevail, the Fourier heat conduction equation for the 
droplet, eq (2), and the boundary condition (9) are replaced 
by a single heat transfer equation. The temperatures in the 
zero-gradient model are denoted by T dz (t) and Tgjjr, t) for 
the droplet and gas respectively. For the case in which the 
droplet temperature is spatially uniform, T d (r, t) = T dz (t) = 
T g (R + , t) = T gz (R + , t) and the heat flux into the droplet, 



(47r/?V3)^c/ 7 f^/(47r/? 2 ) 



(18) 



<£>(77, ct) = /" exp (-0-7)6(77, r)dr 



(26) 



where cr is the dimensionless Laplace transform variable. 
The Bromwich integral 



6(tj, t) 



lim : I 1 

5^cc 2m J y - i8 



exp (ot)<I>(t7, (T)dci (27) 



expresses the temperature in terms of the Laplace transform. 
The quantity y is chosen sufficiently large so that the integral 

J? exp (-<7t)|0(t7, t)|* 

exists. 
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4. Finite-Gradient Model 

From eqs (2)-(9), the Laplace transform of the reduced 
temperature ©(17, r) satisfies the following equations and 
conditions: 



idf) 2 



+ 



d-rf 



5»U - <X 
2 )bf] 



<t> d (r), a) = 



(0 < T) < l) 

(28) 



2 _ orVrA d 
-q 2 /arj 



OrV 



+ o? = 



^(•17, Or) 

(t?>D 



gggtO, Q-) 

3t) 



= 



(29) 



(30) 



[fe + 






[a* 



■a(a»o« 



3 W 

2 Ve£ 



+ 



>M£, o-) = 



(34) 



«*£ 



+ cr- x = 



- ^\*JLc?l, 



o") 



(35) 



where <D d (7}, cr) = ^ d (£, a) and * ff (Tj, cr) = V g (o?t, &). 

The operator in square brackets is the confluent hypergeome- 
tric differential operator and eqs (34) and (35) are confluent 
hypergeometric differential equations. We obtain the general 
solutions to eqs (28) and (29) as 

<J> d (7,, a) = AM (or/6, 3/2, eif/4) (36) 

+ BU(a/e, 3/2, €7, 2 /4) 



for the droplet and 
4> 9 <tj, o-) = o- 1 + 



CW(o/e, 3/2, a*«f/4) .„_. 

+ DU(a/e, 3/2, a*ei»V4) ^ 7) 



lim $(,(17, cr) 

7f-*0O 



4> d (l", cr) = <D B ,1 + , cr) 



d<& d (l- a) a<D/l+ o) 



3i7 



= /3- 



drj 



2o-' 



(31) 



(32) 



(33) 



To solve the subsidiary differential eqs (28) and (29), we 
make a change of variable. Letting £ = er] 2 /^ gives 



for the hot surrounding gas. The quantities M(a, b, z) and 
U(a, b, z) are two independent solutions to the confluent 
hypergeometric differential equation and are called confluent 
hypergeometric functions. Their most useful general proper- 
ties may be found in Chapter 13 of Ref. [11]. 

The constants A, B, C, andZ) must be determined from the 
conditions (30)-(33). The condition that no heat flux exists at 
the center of the droplet requires that B = since U((r/€, 3/ 
2, €T7 2 /4) is not well-behaved at i) = 0. Similarly, the 
condition that the reduced gas temperature goes to unity at 
large distances requires C = since M((r/€, 3/2, a 2 €7) 2 /4>) is 
not well-behaved at infinity. The two remaining conditions at 
the interfacial boundary 7) = 1 constitute two inhomogeneous 
simultaneous equations for the remaining constants A and D. 
Solving these, we obtain the Laplace transformed solutions 
for the reduced temperatures: 



$d(T?, cr) = 
for the droplet, and 



r , id 5 o^e\ {<T 3 aV 



6 2 4 



M4- 



erf 
€2 7 4 



<S> 5 (t7, or) - a" 1 



for th 



e gas where 



o^Hia/e) 



,o- 5 e\ /cr 3 e\ 

2oM[-+ 1, - - + 3XeM - - - 

2 4/ €2 4 



l e 2 4 



o^Hia/e) 



/ / x /°" 5 €\ /cr 3 c^eX , /cr 3 €\ /cr 5 

//(cr/e) = 2M -+ 1, -, - K/ - - — + 3ct 2 /W - -, -)U\-+ 1, - 
7 \€ 2 4/ \€ 2 4 / \e 2 4 / \€ 2 



cr 3 €\_/cr 5 c^e 

4 



(38) 



(39) 



(40) 
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The differential relations of M and U have been used in 
deriving these equations. These expressions have been de- 
rived in a form suitable for going to the limit of zero evapora- 
tion rate constant, € = 0. The correct limiting cases may be 
found in Ref. [6] and the details of the limiting process in 
appendix A. 

The reduced temperature distributions of the droplet and 
surrounding gas are obtained by evaluating the Bromwich 
integrals of eq (8): 



l n 

0(77, r) = lim I y- i8 exp (crT)<&(r), a)dcr. 

1m J y 



8->oo 277"/ J y-i8 



(41) 



To evaluate the Bromwich integrals it is first necessary to 
investigate the nature of the integrand exp (ot)4>(t7, ct). No 
distinction will be made in this discussion between the 



droplet and the gas as the integrands in both cases exhibit the 
same general properties. First of all, exp (ot)4>(t7, ct) is a 
single-valued, analytic function of cr except at its singulari- 
ties. It has a simple second order pole at cr = and an 
infinity of real, negative simple first order poles at the zeroes 
of H{cr/e) given in eq (40). To show that//(cr/e) possesses no 
complex, imaginary or positive zeros is straight-fo reward. 
General methods for locating such zeros may be found, for 
example, in Carslaw and Jaeger [12] and Ince [13]. 

We now evaluate the Bromwich integral. We complete the 
contour of eq (41) by a large semicircular arc in the left-hand 
complex plane and take y to be an arbitrarily small positive 
number. The final reduced temperature distributions are then 
obtained by applying the Cauchy Residue Theorem and 
letting 8 and the radius of the large semicircular arc tend to 
infinity. The integral along the large semicircular arc van- 
ishes and we find 



6,(77, r) 



//(()) 



X dH(p) 
//(O) dp 



- k- 



3 a 2 e 
*2'T 



dp 



dp 



5 a 2 e ', 
+ a 2 pU( 1, -, 



(42) 



+ 3£ 



J . / 5 cre\ / 3 cr 

tp {ep N T)\ p n a 2 pUlp n 4- 1, - , —I - Wipn, -, — 



M\Pn,- ,er) 2 /4 



p\dH(p)/dp\ p=Pn 



e (ri, t) = 1 + 



3x a»(p) 



1 



//(o)L//(o) dp 



3 € 

dMlp, -, - 

H 2 4 



3A- 



p=0 



dp 



-3A- 



BU[ p, - 



3 a 2 er)' 4 



2 4 



dp 



, -, -I 



p=o 



] 



(43) 



I- 



P 2 ndH(p)/dp\ p=Pn 
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where H(p n ) = 0, n = 1, 2, 3 • ••;i.e., the p n are the zeros of 
H(p). 

We pause to make several comments about the expressions 
obtained for GrffTj, t) and 0^(t?, r). Both contain a term 
proportional to the reduced time T which is a consequence of 
the second order pole at cr = 0. As mentioned after eq (9), 
the choice of linear surface regression rate implies that the 
difference between the thermal conductivities times the tem- 
perature gradients evaluated at the interface goes to infinity 
at large times. By large times of course, we mean times close 
to the end of the droplet's lifetime. These are times for which 
the linear surface regression rate of eq (1) is no longer an 
accurate representation of the true physical situation. For A., 
e, and T of the orders of magnitude of interest here, the linear 
term in T may be considered an artifact of the model and 
makes only a very small contribution to the reduced tempera- 
tures. We remark also that eqs (42) and (43) are not suitable 
for going directly to the limit 6 = 0. Further discussion of this 
point may be found in appendix A. 

In section 6 are presented plots of 0^(1, T") versus T for 
typical values of the dimensionless quantities k and e. In the 
next section we derive expressions for the reduced tempera- 
tures for the zero-gradient model. 

5. Zero-Gradient Model 

The analysis for the zero-gradient model is mathematically 
very similar to that for the finite-gradient model and will not 
be quite as detailed. Equations (20)-(25) are the working 
equations for the zero-gradient model. The Laplace trans- 
formed system of equations is 






2 


a 2 er)\ d 
2 )dj) 


<P gz (r), cr) + a 2 


= (44) 




lim <JV(7?, cr) = 1 

17 — >oc 


(45) 




®dz(<r) = <J>« 


(1 + , cr) 


(46) 



for the droplet, and 
<*Wi7, o-) = o- 1 



(2a- + 3\e)<7(cr/e, 3/2, <* 2 eTj 2 /4) 



cr 2 /(o-/e) 



(51) 



for the gas where 



/(o-/e) = 2U(a/e, 3/2, a 2 e/4) (52) 

+ 3/3a 2 c7(l + o/e, 5/2 , a 2 e/4) 

These expressions have been derived in a form suitable for 
going to the limit € = 0, which represents the nonevaporative 
problem solved in Ref. [6], That they in fact reduce to the 
correct limit is shown in appendix A. 

All that remains to be done is to evaluate the Bromwich 
integral of eq (8) with <I>(t7, cr) = O^tj, cr) 

Qgziv, t) = lim — t exp {(n)^ gz {i), (r)dcT. (53) 

The same general comments made about the corresponding 
integrands in the finite-gradient model are applicable here 
also. We complete the contour by a large semicircular arc in 
the left-hand complex plane, apply the Cauchy Residue 
Theorem and let 8 and the radius of the semicircular arc tend 
to infinity. Following this sequence of steps were find 



e d ,(r) = e S2 (i + , t) 



(54) 



and 



QJv, t) = l + 



3 A dl(p) 



1 



/(o)U(o) dp 



a d<J> ,(l + , cr) Xe ,,„, 

- <D dJ (cr) = 3Z \ ~ —■ (47) 

3 dr) 2cr 

The general solution to the differential eq (44) which satisfies 
the condition (45) is 

***(•*?, o) = cr" 1 + AU(<r/e, 3/2, a 2 €7) 2 /4) (48) 

Combining eqs (46) and (47) to eliminate <l> dz (o) yields the 
relation 



ZWi-t&g*. •£ ,49, 

3 or\ 2cr 



which is then used to determine the single unknown constant 
A of eq (48). Doing this yields the Laplace transformed 
reduced temperature distributions 



- 2 - 3X 



dU(p, 3/2, a 2 €7} 2 /4) 



dp 



3eXr 



p = 



(55) 



_ * exp (€ Pn r)(2p n + 3k)U(Pn 3/2, q 2 €T7 2 /4) 
p 2 n d/(p)/dp| p=Pn 



**(<r) = ®M\ o) 



(50) 



17=1 



where I(p n ) = 0, n = 1, 2, 3, • • • ; i.e., the p n are the 
negative real zeros of /(p). Note that the linear term in T also 
appears in the zero-gradient model solutions. 

It is possible to obtain some very simple and accurate 
approximations to eqs (54) and (55). The thermal diffusivity 
of the droplet is several orders of magnitude smaller than that 
of the surrounding gas. This suggests that before applying 
the Laplace inversion formula we expand the subsidiary solu- 
tion (Pgz (77, cr) into a power series in a = aja g . We obtain 
only the simplest approximation by letting a go to zero in eq 
(51): 
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lim <J> 9Z (tj, o-) - <D 90 (i7, a) 

a->0 



(56) 



(T 



2cr + 3eX 
2r)o{cr+ 30) 



The inverse Laplace transform of this remarkably simple 
expression yields the following approximations for the re- 
duced temperatures in the zero-gradient model. 



(57) 



e d0 (T) = (l - -^)(l - exp(- 3/3t)) 
for the droplet, and 

0floO7, t) = 1 - J)' 1 exp(-3j8r) 

- (Xe/2^)(1 - exp(-30$J (58) 



for the surrounding gas. Note that the term linear in r has 
disappeared to the order of approximation used. However, if 
^oziVf °") is expanded out to first order in a, then a rela- 
tively simple infinite series and the linear term in r show up 
in the approximation of Q gz (y, T). Unfortunately, the expres- 
sions obtained by letting a go to zero in the finite-gradient 
model are much the same as the exact solutions. 

In the next section, we present plots of Qdz( T ) versus r for 
typical values of A. and 6. In addition, the times predicted by 
Qdi 7 ), T ), ©dzf 7 ?-* T ) ana< ©do( T ) f° r the droplet to reach a given 
temperature are compared and discussed. 

6. Numerical Examples and Conclusions 

In this section, we give some illustrative numerical exam- 
ples for the predictions made by finite-gradient and zero- 
gradient models. The input data for the calculations are the 
four dimensionless quantities a, /3, X and €. We take a. = 
0.057, j8 - 0.186 and X = 0.01, 0.055, 0.10 and e = 0.10, 
0.55, 1.0, independently, so that there are nine cases. These 
values correspond roughly to those of medium weight fuel oil 
in hot air at atmospheric pressure. Typical values of the 
physical parameters K, d, and C which yield the chosen 
values of a and /3 are given in table 1. In table 2, we give the 
values of L(j corresponding to X = 0.01, 0.055, and 0. 10 and 
k e corresponding to e = 0.10, 0.55, and 1.0. The tempera- 
ture difference T g0 — T d o is taken to be 1666.67 K (corre- 
sponding to T d0 = 300 K and T g0 = 1966.67 K) and 
a d 2 and C d are assigned the values given them in table 1. 

Figures la, lb, and lc compare the numerical predictions 



TABLE 1. Typical values of the thermal conductivity K, density d, specific 

heat C at constant volume, and thermal diffusivity a 2 which yield a = ad/a g = 

0.057 and /3 = K„/K ( , = 0.186 



Property 


Droplet 


Gas 


K(J/cm s °C) 

d(g/cm 3 ) 
C(J/g °C) 
a 2 (cm 2 /s) 


1.45 X 1(T 3 

0.90 

2.39 

6.5 X 10 -4 


2.69 X 10- 4 

1.3 x 10- 3 
1.0 

0.20 



Table 2. The values of the droplet latent heat of vaporization per gram, L ( j, 
and the evaporation rate constant, k (> , corresponding to the values of the 
dimensionless parameters A. = Ld/{C ( i(T g() — T<i„)} and € = k ( ./a,| 2 used. 
The temperature difference T 0o - T,, c is taken to be 1666.67 K and a d 2 and C d are assigned the 

values given in table 1. 



A 


Ui/g) 


e 


k e (vm 2 /s) 


0.01 


39.83 


0.1 


6.5 x 10" 5 


0.055 


219.08 


0.55 


3.575 X \0~ 4 


0.10 


398.33 


1.0 


6.5 X 10" 4 



of the finite-gradient and zero-gradient models for oil droplets 
in air for X equal to 0.01, 0.055, and 0.10, respectively. 
0d(l - , r) and G dz (r), the reduced surface temperatures, 
are computed from eqs (42) and (54). It should be noted that 
due to the presence of the linear term in r, neither 
Grf(l - , r) nor Odz( r ) reach steady state during the droplet 
lifetime. Even if that term were absent, O^tj, t) would never 
become uniform throughout the droplet. We must stress at 
this point that the apparent large difference between the 
temperature profiles for fixed A. but different € is deceiving. 
Although shown on the same graph, the curves are plotted 
versus er/2 which is a function of e, i.e., each curve has its 
own distinct time scale. The results have been shown in this 
fashion for clarity only. If they were all plotted on the same 
scale, real time t for example, the entire group of curves 
would lie very close to one another. As can be seen from the 
figures, G^(l _ , r) and &dz( T ) are quite similar. Their differ- 
ence is greatest when er/2 is small. This is to be expected. 
The thermal gradients inside the droplet are largest near the 
beginning of its lifetime and for such times the approxima- 
tions of the zero-gradient model are more suspect. Although 
not plotted, the simple approximation ( h)(t) to B r /^(T), ob- 
tained in the last section, is an extremely good one as long as 
er/2 is neither too large nor too small. The difference be- 
tween O dz (r) and Brfo( T ) is less than about 15 percent in the 
worst case (X = 0. 1 and e = 1.0) and much smaller in the 
other cases for the r's we have considered. Due to its simple 
analytic form, accuracy, and ease of mathematical computa- 
tion, we suspect that Brfo( T ) is potentially a very useful 
expression for analyses of the initial preheat period of drop- 
lets in hot atmospheres. 

A quantity of interest in the design of combustion cham- 
bers is the time required for fuel droplet to reach a given 
temperature T s . Let us denote this time by t e and consider a 
fuel droplet with the same physical parameters given in table 
1 and an initial radius r^o = 50 jLtm. For illustrative pur- 
poses, let us assume that T$l = 550 K, the initial droplet 
temperature T d o — 300 K, and the initial gas temperature T g0 
= 1966.67 K. The droplet reading T s corresponds to S = 
0.15. Using the data displayed graphically as figs, la, lb, 
and lc, we compute the times required for G to increase from 
zero to 0.15. The results are shown in figure 2 as a plot of t e 
versus e. It can be seen that the times, t e , predicted by the 
finite-gradient and zero-gradient models are within about 20 
percent of each other. For comparison, the t e 's predicted by 
0do( T ) ar e also presented and are seen to be within 10-20 
percent of those predicted by Odz( T ) as long as e < 0.55. 
Figure 2 also provides an extremely powerful cross check on 
the analytic and numerical results of this paper and those of 
Ref. [9]. By solving a totally different problem, Bennett and 
Kayser [9] obtained the times t e for the case of e = 0. These 
are the points on the £-axis of figure 2. By extrapolating the 
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FIGURE 2. The times t e required for a droplet to reach obtained from 
FIGURE 1. The reduced temperature d (7~, t) (solid lines) and B dz (r) d (7~, t) (solid lines and closed circles), O dz (r) (dashed lines and crosses) 
(dashed lines) plotted as a function ofer/2 with a = 0.057, /3 = 0. 186 and e plotted as a function of efor various values of X. 

= 0.10, 0.55, and 1 .0 for X equal tO (a) 0.01, (h) 0.055, and (c) 0.10. The physical properties used are those of table 1 and the initial droplet radius, r d0 , is 50 pun. 

The initial droplet radius, r d0 , is 50 /urn. The temperature T s corresponds to 0{t e ) = 0. 15. 
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results obtained here for nonzero € to e equal to zero, it can 
be seen that the results are in agreement and this indicates 
that this work and that of Ref. [9] are internally consistent. 

Figure 3 contains a plot of the volume fraction of the 
droplet that has evaporated away by the time t — t e versus e. 
For € = 0.1, an evaporation rate constant fairly characteristic 
of the initial preheat stage, the volume fraction evaporated is 
less than 5 percent. If € is fairly small, this perhaps justifies 
neglecting evaporation and surface regression rate in analyz- 
ing the initial droplet preheat stage. A fairly decent approxi- 
mate procedure for predicting the times t e and the droplet size 
at t = t e is to neglect evaporation, compute the time required 
for the droplet to reach a given temperature, and then to use 
that time and a relation like eq (1) to calculate the volume 
fraction evaporated. 
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7. Appendix A 

In this appendix, we show that intermediate results such as 
those of eqs (38) and (51) reduce to the correct limits as e, the 
reduced evaporation rate constant, goes to zero. We also 
make some general remarks concerning this limiting process. 

The results of the 6 = problem corresponding to eqs (38) 
and (51) will be denoted by a superscript '0' and are [6]: 



0.4 



j X= 0.100 



X= 0.100 




FIGURE 3. Volume fraction of a droplet that has evaporated in the time t e 
plotted as a function of efor various values of k. 

The closed circles and solid lines are for the finite-gradient model (6 rf (l~, T)) and the open circles 
and dashed lines are for the zero-gradient model (6*(t)). The physical properties used arc those in 
table 1 and the initial droplet radius, r dih is SO fim. The times t e are those of figure 2. 



•Sfo* «r) 



j8(l + a<r 112 ) sinh(o- ,/2 7)o) 



o-T>o[{j3(l + acr 112 ) - 1} sinh o-" 2 + a 1 ' 2 cosh a" 2 ] 



(Al) 



for the droplet in the finite-gradient model and 



*S*(t?o, a) = o--' - exp[- ao-^7,,,- !)]/[ Vo (* + 3/tacr 1 ' 2 + 3j8)] 



(A2) 



for the gas in the zero-gradient model. The quantities t) and 
r are defined in section 3. We use the relations [11] 



dM(a, 6, z)/dz = aM(a + 1, b + 1 , z)/b 
dU(a, 6, z)/dz - aU(a + 1, 6 + 1, z) 



(A3) 

(A4) 
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to rewrite eqs (38) and (51) as 



®d(ri, or) 



H<-l 



a 3 a 2 er)' 
4 



17 = 1 



(cr 3 a z e\l (cr 3 er)'\ 
\€ 2 4 /J U 2 4 / 



<x 



drj 



/o- 3 eiA /cr 3 a 2 e\ _ /a 3 eW /cr 3 a 2 er? 2 \ 

\7' 2~' T/ _ \ <T' 2~' 4/ ^ \ r'2~ '4 /irT U '2 _ ' ~4 / 



17 = 1 



(A5) 



* w (ij, a) = cr" 1 + - 



(2<r + 3eA)*/(o/e, 3/2, a 2 eri 2 /4) 



2<r[3p(d/dT))U(<r/6, 3/2, a 2 €7j 2 /4)|„ =1 - o-c/(o-/e, 3/2, a 2 €/4)] 



(A6, 



Again from Ref. [11] we obtain 
lim [M(a, 6, z/a)/T(b)] 

a— »oc 

= z 1/2 - 1 ' 2 V 1 (22 1 ' 2 ) 
lim [Hi + a- b)U(a, 6,2/0)] 

a— >oc 

= 2 z " 2 - 1/2 K 6 _ 1 (2z" 2 ) 



(A7) 



(A8) 



where/ and K are modified Bessel functions. Using these two 
limiting expressions, we find 



lim M(cr/€, 3/2, erf/4,)lT(3/2)] 

e+0 



2 sinh(o- 1/2 -n ) 



(A9) 



lim |"lY- - -Wcr/e, 3/2, a 2 €7? 2 /4) 

_ 2V77 r exp(— ai7oO" 1/2 ) 



(A10) 



a^oo- 1 



Multiplying the numerator and denominator of eq (A5) by 

Tf — J/r(— J, substituting the results of eqs. (A9) and 

(A 10), replacing d/dr) by d/drj , and rearranging the result- 
ing expression yields eq (Al); that is 



lim ^(77, cr) = OSXtjo, <*")• 



(All) 



Multiplying numerator and denominator of eq (A6) by 
if —J and then performing the same sequence of opera- 
tions as for the finite-gradient model leads to eq (A2): 



lim <$> 9 



(T7, o-) = ®Wno, °)- 



(A12) 



The same general procedure may easily be used to show that 
other intermediate results also tend to the correct limit as € 
goes to zero. 

It is the intermediate results which may be checked and 
not the final results obtained by evaluating the Bromwich 
integrals. It is not permissible, as may be verified by trying 
it, to take the expressions G^t/, t) and ©02(17, t) of eqs (42) 
and (55) to the limit of e being zero to obtain 
©JK 7 ^ T o) and 6g Z (r)o, t ). This limiting process is im- 
proper for the following reason. Inspection of eqs (Al), (A2), 
(A5) and (A6) shows the first two to be double-valued func- 
tions of cr with a branch points at cr = and the last two to 
be single-valued analytic functions of cr. As one proceeds to 
the limit of zero e, an infinity of simple first order poles 
coalesces into a branch point at the origin, o" = 0. The 
evaluation of the Bromwich integrals leading to ©^(17, t) and 
©02(77, t) is thus only valid for nonzero e. The expressions 
representing 6^(17, t) and O^tj, t) only do so for nonzero e 
and this explains why it is impermissible, in fact impossible, 
to obtain OJ} (i7 , T o) and Qg Z (7) , r ) from them by any lim- 
iting process. 
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